B — ► X s r + r in a two Higgs doublet model 

Yuan-Ben Dai 1 , Chao-Shang Huang 1 , Han- Wen Huang 2 ' 1 
Institute of Theoretical Physics, Academia Sinica, P.O.Box 2735, 
Beijing 100080, P.R.China 
2 CCAST (World Laboratory), Beijing 100080, P.R.China 



Abstract 

The inclusive rate and forward- backward asymmetry of dilepton angular distri- 
bution for a B-meson to decay to strange hadronic final states and a t + t~ pair in 
a two Higgs doublet model are computed. In particular, contributions of neutral 
Higgs bosons to the decay are included. QCD corrections to the effective Hamilto- 
nian for B — > X s t + t~ are calculated using the leading logarithmic approximation. 
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Flavor changing neutral current (FCNC) transitions B — ► X s ^f and B — > X s l + l~ provide 
testing grounds for the standard model (SM) at the loop level and sensitivity to new physics. Rare 
decays B — > X s l + l~(l = e,fx) have been extensively investigated in both SM and the two Higgs 
doublet models (2HDM) [1-13]. In these processes contributions from exchanging neutral Higgs 
bosons can be safely neglected because of smallness of ^-{l = e, //). Recently, r polarization 

asymmetry in B — * X s t + t~ in SM has been studied [jl4| by using the method same as that for 
B — ► X s l + l~(l = e, //). The author of ref.[|l4| did not consider the role played by the neutral Higgs 
boson because the contributions due to exchanging the neutral Higgs boson in SM can also be 
neglected, comparing with those due to 7, Z. 

In this note we investigate the inclusive decay B — * X s t + t~ in a 2HDM. We consider the 
2HDM in which the up-type quarks get masses from Yukawa couplings to the one Higgs doublet H% 
(with the vacuum expectation value V2) and down- type quarks and leptons get masses from Yukawa 
couplings to the another Higgs doublet Hi (with the vacuum expectation value vi). Such a model 
occurs as a natural feature in supersymmetric theories. The Higgs boson couplings to quarks and 
leptons depend on the ratio tg/3 = ^ which is a free parameter in the model. Constraints on tg(3 



from K — K and B — B mixing, T(b — * s^),T(b — * ctu t ) and Rb have been given IE, 16] 



0.7 < tgfi < 0-6(^g) (1) 



(and the lower limit mn± > 200Gev has also been given in the ref. |fL6| ). It is obvious that the 
contributions from exchanging neutral Higgs bosons now is enhanced roughly by a factor of tg 2 (3 
and can compete with those from exchanging 7, Z when tg(3 is large enough. 

Inclusive decay rates of heavy hadrons can be calculated in heavy quark effective theory 
(HQET) [|l7| and it has been shown that the leading terms in 1/mq expansion turns out to be 
the decay of a free (heavy) quark and corrections stem from the order 1/mg [J18|j . In what follows 
we shall calculate the leading term (the 1/mg correction can be easily added if needed). The 
transition rate for b — > st + t~ can be computed in the framework of the QCD corrected effective 
weak hamiltonian, obtained by integrating out the top quark, Higgs bosons and W^, Z bosons 

AC 10 10 

Heff = -^-V tb V t * s (J2 C i^)°i^) +E C Q»W) ( 2 ) 

V 2 i=l i=l 

where Oi(i = 10) is the same as that given in the ref.Q, Qi's come from exchanging the 



neutral Higgs bosons and are denned by 



e 2 



Ql = 16^ (s "™ (fT) 



e 2 



Q2 = ^sffib%)(T7BT) 
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<?4 = yJ^&sxE&gf 
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^5 = ^(^XEfe) 

Q 6 = 7&(^)(E<fc) 

9 

Qio = ^C^^Efe^) ( 3 ) 

with e and 5 being the electromagnetic and strong coupling constants respectively. 

At the renormalization point \x = mw the coefficients Q's in the effective hamiltonian have 
been given in the ref.pl] and Cq.'s are (neglecting the 0(tgf3) term)S. 

C^imiy) = — ~—tg 'p . „ 7 {(aw a + hcos a)fi{x,y) 



m 



2 sin 2 6 w 4 



+ [m 2 /m 2 + {sin 2 a + /icos 2 a)(l - z)]f 2 (x, y) 
sin 2 2a 2 (m^o+m^ ) 2 

+ 2^r (m ^ — 2^ — )/s(y)} ' (4) 



1 We use the Feynman rules given in ref.|l£ 
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where 



CqMw) = -^tg 2 p{h(x,y) + (1 + "^/"^ )h(x, y)}, (5) 
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CQ 3 (m w ) = ^(C Ql (m w ) +C Q2 (m w )), (6) 

o 

CQ 4 (m w ) = —^(Cq^ttiw) ~ C Q2 (m w )), (7) 

C Qi (m w ) = 0, i = 5, ■ ■ ■ , 10 (8) 



m W m H ± 



x , m? 
z = -, h — 



xlnx ylny 

h\x,V) = 7 r, 

x — 1 y — 1 

xlny Inz 
J2{x,y) = t T7 -r + 



(z-x)(x-l) (z-l)(x-ly 

My) = {y _ 1)2 ■ 

Neglecting the strange quark mass, the effective hamiltonian (|2|) leads to the following matrix 
element for b — * st + t~ 

M = ^V tb V t t[Cl ff s LllM b L f^r + C 9 s L ^b L f^ 5 r 

+ 2C 7 m b s L ia tM '^ 1 b R f 1 tl T + C Ql s L b R fT + C Q2 s L b R fj 5 T], (9) 



where 



Ct ff = C 8 + { g (^,s) 



+ 73* 2^ M2 _ n 2_ S M„T, ^ K3C! + C 2 ), (10) 



v ^,M^-q^-iM Vi T Vi 



with s = q /rrifr, q = (p T + + p T -) . In (10) g(^,s) arises from the one-loop matrix element 

of the four-quark operators CVs and has been given in Q. The final term in (10) estimates the 
long-distance contribution from the intermediate tp' fllPf . 
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The QCD corrections to coefficients Cj and CQ t can be incooperated in the standard way by 
using the renormalization group equations. The mixing of the operators Oi(i = 1, 2, • • • , 10) at a s 
order has been studied and the corresponding anomalous dimension matrix has been given j|, |l3| . 
(The mixing of Oi at the next-to-leading order has also been studied in ref. fl3|| .) We studied the 
one-loop mixing of the new set of operators Qi and found that the corrections to coefficients due 
to mixing with Qi are small and can be neglected provided that tg/3 < 50. For example, including 
the mixing between 7 and Q3 (and neglecting the mixing of the other Qi with O7 and the mixing 
of Q3 with the other Qi), we have 



C 7 (m b ) = V ~^{C 7 (m w ) - [^ 10/23 - 1) + ^(v 28/23 ~ l)]C 2 (m w ) - ^ 17/23 ~ WqMw)} 

(11) 

where 77 = a s {m^) / a s {my^) . The additional term is less than 5% correction to the value of C 7 {m b ) 
provided tg/3 < 50. If we include the mixing of Q%(i = 3, 5, 7, 9) with O7 
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(12) 



(and neglect the mixing of Qi(i = 4,6,8,10) with O7) and also consider the mixing among 
Q3) Q5> Q7, Q9 which can be derived from Eq.(39) of ref.p^] and is given by 
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(13) 



with (3q = 11 — |nj, then we have 

58 2Q 

C 7 (m b ) = r ™l™{C 7 (m w ) - [^(r? 10/23 - 1) + Y^(V 28/23 ~ l)]C 2 (m w ) - 0M2C Q:i (m w )} (14) 



The final term in ( |l~4| ) is still small (3% correction) when tg/3 < 50. Using (11) or (|l4[), the limit 



on tg(3 obtained from T(b — > S7) is < 0.25Gev 1 . 
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Qi(i = 1, • • • , 10) does not mix with 0$, Og so that C 8 and Cg remain unchanged and are 
given in ref.||] 



C 8 (m b ) = C 8 (m w ) + ^-^[-^(l- ?? - 11/23 ) + ^(l-r ? - 29 / 23 )]C7 2 (m iy ), (15) 
C 9 (m b ) = Cg(m w ). (16) 



It is obvious that operators Oi(i = 1, • • • , 10) and Qi{i = 3, • • • , 10) do not mix into Q\ 
and Q2 and also there is no mixing between Q\ and Q2- Therefore, the evolution of Cq 1 ,Cq 2 is 
controlled by the anomalous dimensions of respectively. 

C Qi {m h ) = rT^CQ^mw), i = 1,2, (17) 



where 7q = —4 [21] is the anomalous dimension of sl^r. 



Prom Eq.(^), it is easy to derive the double differential distribution d g d J osd as follows 



dsdz y 4vr2/(^) V ' V s \V cb \ 2 

{^ 8 e// | 2 [(l + «)-(l-5)(l-^)z 2 + 4t 2 ] 

1 1 /1+2 

+ eic^tu + ^-a-— )(i- 

s s s s z 

+ ^|C 9 | 2 [(l + s)-(l-5)(l-^)z 2 -4t 2 ] 



2' 1 LX y v /v s 
+ \2Re(Cl ff C* 7 ){l + ^-) + 6#e(C 8 C 9 *)(l - ^-)hz 



Af2 , Q Q 

+ URe{C 7 C* 9 ){l - —)*z + -\C Ql \\s - 4i 2 ) + -\Cq 2 \H 
At 2 , 4t 2 - 



+ 6i?e(C 8 C^)(l - — )*tz + l2Re{C 7 C* Ql )(l - —)*tz 
+ 6Re(C 9 C* Q2 )t}, (18) 

where z = cos9, t = — and 6 is the angle between the momentum of the B-meson and that of 
/ + in the center of mass frame of the dileptons t + t~ . Integrating the angle variable, we obtain the 



5 



invariant dilepton mass distribution 



dT(B -> X s t+t-) a 2 4t 2 g |y# t ;p 

_ ^^x>) i -^ y (i- s) (i--) 2 ^ F - 

{\C! ff \ 2 (l + ^)(l + 25) + 4|C 7 | 2 (1 + ^)(1 + f ) 

2+2 2/ 2 
+ |C 9 | 2 [(1 + 2s) + —(1 - 4s)] + 12Re(C 7 Cl* ff )(l + -j-) 

+ l\C Ql \ 2 (s-4t 2 ) + ^\C Q2 \ 2 s + 6Re(C 9 C* Q2 )t}. (19) 



Switching off the channel exchanging neutral Higgs bosons, Eq,(|l8|) turns out to be the same as 
that given in Ref.[14]. We also give the forward-backward asymmetry in B — > X s t + t~ 

f\ d2 v f° J d 2 r 

A{s) = / dz— — - / dz— — 

Jo dsdz J-x dsdz 



{Re(C 8 C;)s + 2Re{C 7 Cl) + Re(C s C* Ql )t + 2Re{C 7 C* Ql )t] . (20) 

The following parameters have been used in the numerical calculations: 

m t = 175Get>, = 5.0Gev, m c = 1.6Gev, m T = 1.7Gev, A-p^ = O.lGev. 

Numerical results are summarized in Fig. 1-2. 

It is obvious that there are sharp peaks on the curves of differential branching ratio and 
asymmetry A(s) around the mass of resonance ip'. One can see from the figures, that the contribu- 
tions of neutral Higgs bosons is significant in 2HDM when tg/3 is larger than 25 and the values of 
the masses of the Higgs bosons are in the reasonable range, and the forward-backward asymmetry 
of dilepton angular distribution is more sensitive to tg/3 than the invariant mass distribution. We 
expect that the more precise experiments on rare decays of B mesons will shed light on the two 
Higgs doublet model. 

This work is supported in part by the National Natural Science Fundation of China and 
the Postdoctoral Science Foundation of China. 
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Figure 1: The differential branching ratio (normalized to T(B — > X c eu e )) as a function of the 
scaled invariant dilepton mass s = s/mf in the decay B — > X s t + t~. We have taken m^± = 200Get>, 
m^o = 80Gev, m#o = 150Get>, and m^o = lOOGew. Assumed tg/3 = 30,50 are indicated on the 
curves. The solid line corresponds to the case of switching off the channel exchanging the neutral 
Higgs bosons. 



— I 1 1 

tg(3 = 30, C Qi = 

tg(3 = 30 

tgP = 50 ■ 




Figure 2: Forward-backward asymmetry of the dileptons on the decay B — > X s t + t~ , A(s) (nor- 
malized to T(B — ► A" c ez> e )), as a function of the scaled invariant dilepton mass s = s/mf. We have 
taken m H ± = 200Get>, m ft o = 80Gew, m#o = 150Gev, and m A o = WOGev. Assumed tg/3 = 30, 50 
are indicated on the curves. The solid line corresponds to the case of switching off the channel 
exchanging the neutral Higgs bosons. 
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